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ABSTRACT We propose a new method for carrying out field-theoretic 
simulations of polymer systems under conditions of prescribed external stress, 
allowing for shape changes in the simulation box. A compact expression for 
the deviatoric stress tensor is derived in terms of the chain propagator, and 
is used to monitor changes in the box shape according to a simple relaxation 
scheme. The method allows fully relaxed, stress free configurations to be ob- 
tained even in non trivial morphologies, and enables the study of morphology 
transitions induced by external stresses. 

last modified February 2, 2008 

1 Introduction 

In particle-based simulations, such as Molecular Dynamics (MD) or Monte 
Carlo (MC), it has long been realized that the flexibility of the approach 
could be greatly improved by imposing intensive, rather than extensive, 
thermodynamic constraints. Schemes to impose pressure are well known in 
Monte Carlo simulations, and constant temperature and constant pressure 
algorithms, using various thermostats and barostats, have become common 
options in Molecular Dynamics simulations £Q. When dealing with solids, 
however, it was pointed out by Parrinello and Rahman [2] that the study 
of structural transformations would necessitate a release of the constraints 
imposed on the symmetry of the simulation cell. They introduced, within the 
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context of MD simulations, a method that treated the cell shape as a dynami- 
cal variable, coupled to an externally imposed stress tensor. The method was 
later slightly modified by Ray and Rahman [3| to account for some subtleties 
associated with describing large deformations. The Parrinello-Rahman-Ray 
technique has been used in a variety of contexts, e.g. to describe structural 
transformations between crystalline phases in solids and to compute elas- 
tic constants under pressure, and is still evolving as new applications are 
identified E] • 

In the polymer literature, such variable cell shape methods have been ap- 
plied in particle based simulations (see e.g. [7j). However, in the other promi- 
nent class of simulation methods for complex fluids, i.e. the "field-theoretic" 
approaches [Sj, no equivalent of the variable shape simulation cell methods 
has been proposed up to now. The most familiar type of field-theoretic sim- 
ulation is a numerical implementation of self consistent field theory (SCFT) 
[HI E3> which amounts to the computation of a saddle point (mean-field) 
configuration of a statistical field theory model. Modern SCFT simulations 
are normally of two types. The so-called large cell calculations are used to 
address bulk fluid systems with no long-range order, or microphase separated 
polymeric fluids where the symmetry of the periodic domain structure is not 
known in advance. Cubic simulation cells with periodic boundary conditions 
are normally applied for this purpose [TTJ O E2] and large cells must be 
used to avoid finite-size effects. Often the cell size is adjusted in such cal- 
culations to minimize the the free energy, although the cell shape is fixed. 
In the second type of SCFT simulation, a so-called unit cell calculation, one 
assumes a particular crystallographic symmetry and applies a simulation cell 
that corresponds to a primitive cell compatible with that symmetry [TUllTl] . 
In these calculations, the cell size and shape are adjusted to minimize the 
free energy of a unit cell of the periodic structure No analogous proce- 
dure currently exists for adjusting the cell shape in large cell calculations in 
order to achieve stress-free conditions. 

In view of the growing interest in the computation of the mechanical 
properties of the complex mesophases that can be formed in copolymer melts 
[TEl ll7|. it is of interest to formulate a method similar to the Parrinello- Ray- 
Rahman approach within the framework of field-based polymer simulations. 
Such a "variable cell shape" approach should provide a useful option in SCFT 
calculations when dealing with systems subjected to anisotropic stresses (e.g. 
unaxial extension or shear), or in situations where the optimal symmetry of 
the simulation box is not a priori obvious. Even more simply, it would 
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provide a simple means of relaxing the box size and shape to an optimal 
value, simply by imposing a 'zero stress' external condition. In the present 
manuscript, we present a first step towards such a simulation methodology in 
which the cell shape is adjusted at fixed cell volume. Such a formalism is par- 
ticularly appropriate for field theory models that invoke the incompressible 
liquid assumption. 

The paper is organized as follows. In Section |2l we recall some useful 
notations and results from elasticity theory. The expression for the stress 
tensor appropriate for polymeric systems is discussed in Section EJ where we 
also describe the equations used for evolving the shape of the simulation cell. 
Some examples of two-dimensional simulations for AB diblock copolymer 
melts are given in Section 



2 General Notations 

In this section, we closely follow the description made by Ray and Rahman 
[3] of the Parrinello- Rahman type methods used in particle simulations. In a 
system with periodic boundary conditions, the simulation cell can be defined 
by three (in general nonorthogonal) independent vectors h!,h 2 ,h 3 forming the 
sides of the parallelepiped cell. The Cartesian coordinates of these vectors 
can be used to construct a 3 x 3 matrix h defined by h = (hi, h 2 , h 3 ). The 
Cartesian coordinates of any point R in the cell can be expressed as 

R=hX (1) 

where X is a rescaled vector whose components lie in [0, 1]. Integrals on R 
can be converted into integrals over X by using a scaling factor det h, which 
represents the volume of the cell, V. In the case of a particle or monomer 
number density, for example, one can write 

p(R) = p(X)(deth)- 1 (2) 

The metric tensor G is constructed from h as 

G = h T h (3) 

where is the transpose of h. G is used in transforming dot products from 
the original Cartesian to rescaled coordinates, according to 

R • R' = X • G • X' = X a G a pXp (4) 
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where here and in the following summation over repeated indexes is implicit. 

Elasticity theory describes the deformation of any configuration from a 
reference configuration in terms of a strain tensor. This tensor is constructed 
by relating the vector connecting two points in the deformed configuration 
to the corresponding displacement of the same points in the reference con- 
figuration. If the reference configuration of the simulation box is denoted by 
h , the strain is given by 

e = \ [(h T rh T h(h r -i} = \ [(h T rG(h r - 1] (5) 

where 1 denotes the unit tensor. It is important to note that this expression, 
usually known as the "Lagrangian stress tensor" is not limited to small de- 
formations [TH]. Usually, the reference configuration will be defined as a 
state of the system under zero applied external stress. 

The thermodynamic variable conjugate to this strain tensor, in the sense 
that the elementary work done on the system can be written in the form 

SW = V Tt(t e), (6) 

is the thermodynamic tension tensor r [T^], also known as Piola-Kirchhoff 
second stress tensor. V = det h denotes the volume of the system in the 
reference configuration. This thermodynamic tension tensor can be related 
to the more usual Cauchy stress tensor q_ through 

*=yh(h r Z (h\rh T ( 7) 

The tension is the derivative of the free energy with respect to the strain, 
which is calculated from the reference configuration. The stress, on the other 
hand, is the derivative of the free energy with respect to an incremental 
strain taken with respect to the actual configuration. The difference between 
these two quantities can be understood, qualitatively, from the fact that the 
strain is not an additive quantity, as can be seen from the existence of the 
nonlinear term in eqEl While the Cauchy stress has a mechanical meaning in 
terms of forces applied to the sample, the thermodynamic tension is a purely 
thermodynamic quantity, and does not in general have a simple mechanical 
interpretation. 

It was proposed by Ray and Rahman [3] that, in a particle system, an 
appropriate ensemble to carry out simulations allowing for shape transfor- 
mations in the simulation cell is a constant tension, constant temperature, 
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constant number of particles ensemble, usually denoted as a "rTN" ensem- 
ble. In this ensemble, fluctuations in the shape and size of the simulation 
cell, described by h, are allowed, and the configurational partition function 
is written as 

Z = J dhj d N Rexp(-pV r :e-f3H(R)) (8) 

where R denotes the coordinates of the N particles, which evolve in the 
"box" defined by h, [3 = 1/(/cbT), and H(R) is the potential energy In eq 
IH1 no kinetic energy terms in the Hamiltonian associated with either particle 
or cell motion were included. Such terms are integral in the Parrinello- 
Rahman method in order to obtain equations of motion for the cell and the 
particles. In SCFT theory, however, the inertial particle dynamics is replaced 
by overdamped relaxational dynamics for the density fields, and we will make 
use of a similar relaxation scheme for the cell coordinates. Thus, the kinetic 
energy terms prove to be irrelevant. 



3 Application to Polymer Field Theory 
3.1 Polymer Partition Function 

The key step in formulating a polymer field theory is to transform the integral 
over the particle coordinates appearing in the partition function for a complex 
fluid (c.f. equation|SJ) into an integral over density and/or chemical potential 
fields for the different species. This step has been described in many papers 
(see e.g. [HI EI); so we will only highlight here the technical modifications 
associated with a variable cell shape simulation method. To illustrate the 
method, it is convenient to work with a standard mesoscopic model of an 
incompressible AB diblock copolymer melt ^Uj- In this model, the melt 
consists of n identical diblock copolymer chains composed of monomer species 
A and B and contained in a volume V. Each of the chains has a total of N 
statistical segments; a fraction / of these segments are type A and constitute 
the A block of each macromolecule. For simplicity, the volume occupied by 
each segment, v , and the statistical segment length, b, are assumed to be 
the same for the A and B type segments. The Hamiltonian for this system 
can be written 

^bT f 1 f dHi(s)\ 2 f ^ „ 

H = X^4jp- J ds y dg J + v oXabK B T J drp A (r)p B (r) (9) 
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where Rj(s) with s e [0, 1] is a space curve describing the conformation of 
the ith copolymer and R 2 gQ = b 2 N/6 is the radius of gyration of an ideal 
chain of N statistical segments. Interactions between dissimilar segments A 
and B are described by the Flory parameter xab- The densities Pa,b(j) are 
microscopic segment density fields defined by 



p A (r)=iVV I ds 5(r-R, 
i=i Jo 



») (io) 



and 



n „\ 

p B (r) = iVV / ds5{r-B. t {s)) (11) 
i=i ■ / / 



A local incompressibility constraint Pa( y ) + Pfi( r ) = Po is imposed in this 
standard copolymer melt model for all points r in the simulation domain. The 
total segment density po can evidently be expressed as po = nN/V = l/v . 

In order to implement the Parrinello-Ray-Rahman method, it is conve- 
nient to switch from the Cartesian coordinates R to the rescaled coordinates 
X, as defined in eq^ Using the rescaled coordinates, the Hamiltonian of eq 
El is expressed as 

~t AR 9 o Jo ds ds 

+ v XAB(det h)' 1 J c/xp j4 (x)p B (x) (12) 

where the coordinates X and x are taken in [0, l] 3 . The partition function is 

Z = J d(h)(deth) nN 5(deth- V ) exp(-(3V T : e) 

x J] / VXiew(-PH)]l6tf A W+p B (x)-nN) (13) 

The final factor in the above expression imposes the constraint of local in- 
compressibility for the scaled density fields. Moreover, the assumption of 
incompressibility implies globally that the cell volume remains fixed at its 
initial value, i.e. det/i = V = Vq = deth . This is enforced by the delta 
function in the first line above. All shape transformations considered in the 
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present work are therefore volume preserving. The practical implementation 
of this constraint will be discussed below. 

Hubbard-Stratonovich transformations are used to convert the particle- 
based partition function into a statistical field theory [S] . These can be 
carried out straightforwardly on the polymer partition function for a given 
cell shape h, Z(h), with the result 

n „ 

Z(h) = J] VX t exp(-f3H)l[5(p A (x)+p B (x)-nN) 

i=l J x 

= J Vw exp(nlaQ[w,h]- E[w]) (14) 

where Q[w,h\ is the partition function of a single copolymer chain experi- 
encing a chemical potential field w(x, s), J Vw denotes a functional integral 
over the field w, and E[w] is a local quadratic functional of w that reflects 
the A-B monomer interactions and the local incompressibility constraint: jH] 



dx 



(iu B - ioa) - (w A + IU B ) 



(15) 



_2 X iV 

Here we have noted that for an AB diblock copolymer melt, the potential 
w(x,s) amounts to a two-component potential, i.e. w(pc,s) = iu^(x) for 
s G [0, /] and io(x, s) = to B (x) for s G [/, 1]. 

The object Q[w,h\ is a normalized partition function for a single copoly- 
mer experiencing a potential field u>(x, s) and is defined as 

f DX e ~ 



JQ as — as 



The path integral in the denominator ensures the normalization Q[0, h\ = 1, 
but it depends only on deth = V, so it remains constant for volume- 
preserving changes in cell shape h. The partition function Q can be obtained 
from a single-chain propagator g(x, s) that is the solution of a modified dif- 
fusion equation 



^o(G-%^r4r- - s )<?(x, s ) ( 17 ) 



subject to q(x, 0) = 1. The single chain partition function is given by 
Q[w,h\ = J* dx g(x, 1). 
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By virtue of the above equations, the partition function for an incom- 
pressible diblock copolymer melt confined to a cell of variable shape can be 
expressed clS 8b field theory in the variables h and w. 

Z = J d(h) J Vw 5(deth-V ) exp(-F[w,h]) (18) 

where F[w,h\ is an effective Hamiltonian given by 

F[w, h] = pV r : e + E[w] - n In Q[w, h] (19) 

The factor (deth) nN = Vq , which was present in eq El has no thermo- 
dynamic significance for a constant volume system, so it has been absorbed 
into the integration measure in eq^J This equation provides a rather gen- 
eral starting point for carrying out variable cell shape simulations of polymer 
field theories. Here we shall invoke the mean-field approximation (SCFT) 
and, for a given shape h of the simulation box, approximate the functional 
integral over w in eq^Jby the saddle point method. For this purpose, the 
functional Q[w,h\ can be evaluated for any w and h by solving the modified 
diffusion equation using a parallel implementation of a pseudospectral oper- 
ator splitting scheme [201 El- The saddle point (mean-field) value of w, w*, 
is obtained by applying a relaxation algorithm [HI El E] to solve 



5F[w,h] 



5w(x, s) 



(20) 



w=w* 



In the mean-field approximation, F[w*,h\ corresponds to the free energy of 
the copolymer melt (in units of ksT). 



3.2 Equation of Motion for the Cell 

In the Parrinello- Rahman method, the cell is ascribed a fictitious kinetic 
energy, and the corresponding Hamilton equations are solved along with the 
equations of motion for the particles. As the equations of motion for the 
fields in SCFT are overdamped relaxation equations, a reasonable choice is 
to use a similar relaxational dynamics to evolve the cell coordinates. We 
have implemented the following dynamics for the evolution of h: 

*._^«fc4l (21) 
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where the tensor D_ is a projection operator whose action on an arbitrary 
tensor M is a traceless tensor, i.e. D M = M — (l/3)Tr(M)l. Equation 
|2*T1 corresponds to a cell shape relaxation that (for Ao > 0) is down the 
gradient dF/dh, approaching a local minimum of the mean-field free energy 
F[w*,h\. The "mobility" tensor hPh ~ l is chosen so that the cell shape 
dynamics described by eq IIH] conserves the cell volume. This can be seen 
from the identity 

d , , . , „ _ 1 d 



— In det /i = Tr(h~ l —h) 
dt ~ K ~ dt~' 



(22) 



which implies that det h is a constant of the motion if hT 1 (d/ dt)h is traceless. 

Application of eq |^ requires an expression for the thermodynamic force 
dF/dh. Explicit differentiation, noting the constraint of constant det h, leads 
to 

dF[w,h] _ / d 



dh 



(23) 



where S is a symmetric tensor defined by 



E af3 [w,h] 



2k B Tnd\nQ[w,h] 



V 



dG 



ds 



dX a (s) dXp(s 



ds 



ds 



(24) 



The angular brackets in the second expression denote an average over all 
conformations X(s) of a single copolymer chain that is subject to a prescribed 
chemical potential field w and fixed cell shape h. 

The first term on the right hand side of eq E21 can be conveniently rewrit- 
ten as 







.(Tr(e_r )) 







dh a0 dh a p K '2 

Hence, eq ED can be compactly expressed as 

dh 
~dt 



(Tr-^ T h T h h^ l r ) = (hh-'r h T ) af3 . (25) 



-XhD 



(26) 



where A > is a new relaxation parameter defined by A = /3V()\q. 

Equation |2H1 is a purely relaxational equation that will evolve the cell 
shape to a configuration of minimum free energy (in the mean-field approx- 
imation). This configuration can either be metastable (local minimum) or 
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stable (global minimum). Addition of a noise source to the equation provides 
a means for overcoming free energy barriers between metastable and stable 
states, i.e. a simple simulated annealing procedure. Beyond mean-field the- 
ory, in order to fully assess the effects of thermal fluctuations in both w and 
h, a more sophisticated procedure is required that takes into account the 
non-positive definite nature of the polymer field theory [S]. 

An equilibrium solution of the cell shape equation|2E]is evidently obtained 
when 

(^I^oVs =0 (27) 

Combining eqs El and EH it is easily seen that this equilibrium condition 
corresponds to a balance between the externally applied Cauchy stress, g_ , 
and the internal elastic stress, a mt , sustained by the polymer chains 

a + a M = (28) 

where 

^^meu^t{j\ d ^). ( 29) 

This expression for the internal polymer stress is well-known in the polymer 
literature [22] • 

Equation EH1 drives a change in the shape of the simulation cell (at con- 
stant cell volume) to approach the equilibrium condition EH at which the in- 
ternal elastic stress of the copolymers balances the imposed external stress. 
In practice, our numerical strategy consists in implementing a cell shape 
evolution at regular intervals, typically every 50 to 100 field minimization 
steps. The magnitude of the coefficient A in eq I2H1 should ensure that the 
strain associated with such evolutions remains moderate. As the driving force 
itk) l L h^o ) + ^ 1S normally of order unity (in the appropriate reduced units 
used in SCFT calculations, i.e. taking ksT as the unit of energy and R g0 as 
the unit of length), 1/A is a number of the order of the interval between cell 
shape evolutions. 

A final technical point is that the change in cell shape should ideally 
be performed at a fixed spatial resolution for the solution of the diffusion 
equation. Hence, a change in the cell shape generally involves re-meshing 
the rescaled simulation cell with a grid that will be finer in the expanded 
directions and coarser in the compressed directions, in order to keep the 
mesh size locally invariant. 
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3.3 Expression for the Stress Tensor 

The last step in the implementation is to find an expression for the internal 
stress tensor a mt (eg 129(1 or X f eg 124}) in terms of the single chain propagator, 
which is the central object computed in a field-theoretic simulation [51 12~5]. 
In order to obtain such an expression, we return to the definition of E 

m/v 9 1 dQ[w,G] 

f3VY, a p = -2n^-^— (30) 

The derivative of the single chain partition function can be calculated by 
discretizing the paths with a small contour step A. 

J VX(s)5(X - X(s))5(X> - X( S + A)) (fc^) (^) 
exp Jo' ds d -^G a ^ - £ ds w(X(s), s)) (31) 



Except between the points X, s and X', s + A one can replace the path 
integrals with propagators q, so that 

- = 4^ £ ds Id * 1 dX ' 9(x ' s)g(x ' 1 - s - A) 

(^A^) ( "a S ) eX P ("lA^^W 3 ^ ~ X 'a)( X P _ X^ 

One can then set X' = X + u, and expand for small u and A according 

to 

da da d 2 a 
5 (X + »• 1 — A) = ?(X, 1 - .) - + + (33) 

The derivative w.r.t. s can be eliminated by applying the modified diffu- 
sion eq El One also requires second and fourth moments of the Gaussian 
distribution of displacements u, 

u a u^us = (2R%A) 2 (G- 1 p G^ + G^G^ + G^G^) 

11 



(32) 



By means of these results, we have 

+ ( G ^ G 7 * + G ai G 05 + GjGfa) / rfx Jo dsq(X, s) a ®J xg 

(34) 

where p(X) = Q" 1 Jq 1 c/sg(X, s)g(X, 1 — s) is the single-chain total monomer 
density operator. There is a partial cancellation in the last two terms so that 

- = G~l 12$* 4- iG- / <OC»(X)p(X) 

(35) 

The internal polymer stress is obtained after matrix multiplication by h 
on the left and /i T on the right. This implies that the first two terms become 
a simple isotropic stress contribution, and are therefore not relevant to an 
incompressible system. The final formula for the internal stress tensor is 
therefore, apart from this diagonal contribution, 

Ant (nk B T\ 2R% f /- 1 &q(X,l-s) x 

The tensor £ appearing in equation |2H1 is given by an expression similar to 
I3TH with G replacing h. The factor ksTn/V accounts for the total number 
of chains, and produces a stress with the correct dimensions. In practice, 
the stress will be made dimensionless by dividing by this factor, so that the 
dimensionless stress is given by 



(„ w q « J dX I ds « (x ' " im, h » (37) 

An equivalent expression that better reveals the symmetry of the stress tensor 
is obtained by an integration by parts: 

"3 - ! \f ffl f^». (38, 



{n/V)k B T Q 7Q J J 3X 1 dX, 
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A local (rather than volume averaged) version of this connection between 
the stress tensor and the polymer propagator was derived previously in j2H] . 
Numerically, a a f3 is evaluated from eqEZI using the pseudo-spectral scheme of 
PUI I2T] . The derivatives with respect to spatial coordinates are obtained by 
multiplying the propagator by the appropriate components of the wavevector 
in Fourier space, and transforming back into real space. 

4 Examples 

4.1 Stress-Strain Relation for Lamellae 

As a first example and check of the method, we consider the stress-strain 
relations for a perfectly aligned lamellar system, with lamellae perpendicular 
to the z direction. Such a system can support a stress only in the z direction 
and is therefore described by a single elastic relation, a zz versus e zz (Other 
"apparent" elastic constants computed in J7j can be shown to vanish, as is 
obvious from eq ESI m a system that is translationally invariant along the x 
and y directions). For such a system, we have computed the stress-strain 
relation using either an imposed strain method or an imposed stress, as de- 
scribed above. The imposed strain amounts to progressively deforming the 
simulation cell along the z direction, keeping the topology (i.e. the number of 
lamellae) fixed. The stress is then evaluated using eqEU The reference con- 
figuration was created by running the variable shape simulation code under 
zero stress conditions (zero tension), in order to obtain a perfectly relaxed 
initial configuration. The results, shown in Figure ^ indicate that the two 
methods produce identical results in this simple situation. Finally, we should 
mention that lamella under dilational stress should give rise to a "chevron" 
instability at rather small strains, as described in [21]. In order to avoid the 
development of this instability, the equations were solved without any added 
noise, which allows one to preserve the metastable lamellar topology even at 
large strains. 

4.2 Formation of Cylindrical Phase 

As a second example, we show how a square cell, quenched into the region 
of the phase diagram where the equilibrium phase is a hexagonally-packed 
array of cylinders, naturally evolves under zero tension to a rhombic shape 
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characteristic of a primitive cell of the Bravais lattice. Starting from a random 
initial field configuration and imposing a quench to \N = 15.9 (for a diblock 
copolymer with / = 0.64), successive stages of this transformation are shown 
in Figure El A cylindrical nucleus appears after a few iterations, and the 
cell progressively deforms from a square to a rhombus consistent with the 
triangular organisation of the lattice. Relaxation to this equilibrium shape 
proceeds in typically 1000 iterations, with a change in cell shape implemented 
every 10 iterations. 

4.3 Shift of the Lamellae- Cylinder Transition Under 
Stress 

As a final application example, we attempt to address the following question: 
how does a compressive stress affect the lamellar to hexagonal transition in 
a two dimensional system? We first determine the location of this transition 
for a simple AB diblock copolymer melt with the fraction / of A monomers 
equal to / = 0.64 by running a "zero stress" simulation for various values 
of xN- As the lamellar spacing or triangular lattice spacing change with 
X-/V, the zero stress simulation represents a convenient way of automatically 
adjusting the cell dimensions to the minimum free energy, thus obtaining 
fully relaxed, stress free configurations at each value of %iV. The resulting 
transition takes place at xiV ~ 18, as shown in Figure HU which was obtained 
by ramping xiV at zero tension. 

As the lamellar spacing at the transition is smaller than the spacing be- 
tween cylinders of the triangular lattice, it can be expected that the applica- 
tion of a compressive stress will favor the lamellar phase over the hexagonal 
one. In order to verify this hypothesis, we have carried out simulations in 
which the hexagonal phase at \N = 16 is subjected to an external uniaxial 
compressive stress of 1 in reduced units. Under this imposed compression, 
the initially almost rectangular cell shape transforms into a parallelogram, 
until suddenly the structural transformation takes place (see Figure|3]). Com- 
paring the free energies under tension, it is found that the lamellar phase has 
indeed a lower free energy when the transition takes place. However, since the 
resulting lamella are not aligned with the principal directions of the stress, 
they cannot sustain the resulting shear so that the cell eventually elongates 
without bound. 
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5 Conclusion and perspectives 



Variable shape simulations have played, and are still playing, an important 
role in molecular simulations of solids for the identification of allotropic phase 
transformations under pressure and of the associated transformation path- 
ways. Here we have shown that the methods used in solids can be generalized 
to simulations of polymer field theories. In this context, many types of poly- 
meric fluids can be treated as incompressible fluid mixtures. As a first step, 
we have therefore focused our attention on the incompressible case and de- 
veloped a methodology by which a simulation cell can be adjusted in shape 
to attain the lowest free energy in the mean-field (SCFT) approximation. In 
the case of zero applied stress, the method provides a very convenient way 
of obtaining mechanical equilibrium in large systems, which may include de- 
fects or complex structures. We have also shown that the application of a 
non-hydrostatic tension can modify the thermodynamic equilibrium among 
phases, and that a variable shape simulation cell allows the system to explore 
non-trivial pathways for structural transformations. 

The relaxational sampling of cell shapes that we have used in this paper is 
rudimentary, and there is certainly room for improvement if one is interested 
in studying, e.g., nucleation barriers between phases. Another interesting 
perspective is that, although we have been treating the internal stress ten- 
sor as a global property of the system, it is clear from eqEHl (see also [22]) 
that a local value of the stress tensor can be easily identified prior to volume 
integration. This formula should be useful in assessing the local stress dis- 
tribution in inhomogeneous systems under strain, including polymer-based 
composites. 
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Figure 1: Stress-strain relation for perfectly aligned lamellae, where the stress 
is expressed in units of (n/V^/cgT according to eqEZI The full line is obtained 
from imposed strain simulations, while the dots are obtained using constant 
tension simulations. The dashed line corresponds to the linear elastic behav- 
ior around the zero stress configuration. The associated elastic constant is 
B = 4(p/N)k B T 
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Figure 2: Transformation of a square cell under zero external stress, when the 
melt is quenched into the stability region of the cylindrical phase (\N = 15.9, 
/ = 0.64). 
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Figure 3: Free energies (per chain) of the lamellar (solid) and cylindrical 
(dashed) phases under zero tension, as a function of \N , for a diblock with 
/ = 0.64. 



18 



I 



14 
12 
10 



« 



14 



12 



10 



1 2 3 4 5 



1 2 3 4 5 



14 



• * ' 

HP 

'lit 

r 

1 2 3 4 5 



14 



PPM 

1 2 3 4 5 



! I 

■ w 



14 



10 




14 



12 



10 



1 2 3 4 5 



1 2 3 4 5 



Figure 4: Six stages of the cylindrical-lamellar transformation upon com- 
pression of a cylindrical phase at \N = 16 along the horizontal axis 
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